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Abstract 

We analyze the effect of variation of fundamental couplings and mass scales 
on primordial nucleosynthesis in a systematic way. The first step estab- 
lishes the response of primordial element abundances to the variation of a 
large number of nuclear physics parameters, including nuclear binding en- 
ergies. We find a strong influence of the n—p mass difference (for the "'He 
abundance), of the nucleon mass (for deuterium) and of ^ = 3, 4, 7 binding 
energies (for ^He, ^Li and ^Li). A second step relates the nuclear param- 
eters to the parameters of the Standard Model of particle physics. The 
deuterium, and, above all, ^Li abundances depend strongly on the average 
light quark mass rh = {mu + md)/2. We calculate the behaviour of abun- 
dances when variations of fundamental parameters obey relations arising 
from grand unification. We also discuss the possibility of a substantial 
shift in the lithium abundance while the deuterium and ^He abundances 
are only weakly affected. 

1 Introduction 

The constancy over space and time of the coupling strengths and particle masses 
in the Standard Model of particle physics is an assumption that should be 
testedQ If a variation existed, it would violate the condition of Local Position 
Invariance contained in the Einstein equivalence principle of General Relativ- 
ity: non-gravitational experiments measuring the same quantities at different 
times would give different results. Within a relativistically covariant setting, 
a variation arises due to the coupling of Standard Model particles to a scalar 
field whose cosmological value depends on time [21 13] . 

Dynamical dark energy or quintessence [H [5] is precisely due to a scalar field 
whose value continues to vary in recent cosmological epochs. Thus the possibil- 
ity of time-varying couplings has been discussed since the first suggestions 
of dynamical dark energy. Related ideas concern the dilaton or time-varying 
moduli fields in string theory Many dark energy candidates in higher- 
dimensional or string theories can be excluded because the time-variation of 
couplings is too strong. On the other hand, a detection of time-varying cou- 
plings would be a strong hint in favour of a dynamical dark energy, explaining 
the renewed interest in this subject P [9l [TOl [TTl [T2| [T3] . 

^For a general review of "varying constants" see [I]. 



Such investigations have recently been further motivated by possible signals 
of a nonzero variation at redshifts about 0.5-4, due to astrophysical absorption 
spectra which deviate from those found in the laboratory. By choosing the 
transitions to be measured, one can achieve sensitivity at the 10~^ level to 
the fine structure constant a, the proton-electron mass ratio fi = nip/me |14j . 
the proton gyromagnetic ratio gp, and various combinations of these [15]. The 
current observational situation is contradictory, with both positive |16j and 
null [T71 [TBJ results, and continuing discussions of methods and of statistical 
and systematic errors |19j . 

Other limits on the values of particle couplings and masses at (relatively) 
recent epochs arise from nuclear physics effects in the Oklo natural reactor at 
z ~ 0.2 \20\ I21| and from long-lived /? decay isotopes in meteorites |22j. Direct 
comparisons of atomic clocks at periods of a few years in the laboratory have 
also led to strong bounds at the level of 10~^^ per year fractional variation [23]. 
Since these are Earth- or Solar System-based probes their results may not be 
directly comparable with those from cosmologically distant absorption systems, 
without further theoretical assumptions. Clearly it is desirable to have as many 
bounds as possible at different redshifts and in different environments to probe 
possible variations in a model-independent way. 

In this paper we consider possible variations at a much higher redshift, 
that of Big Bang nucleosynthesis (BBN) at z ~ lO^'', which is currently the 
earliest time at which theories of nuclear and particle physics can be compared 
to cosmological observation in a controllable way \24\ [25| [26] . It is remarkable 
that primordial abundances are influenced by every known force of interaction: 
gravity, through the expansion rate of the Universe; weak interactions, through 
neutron decay; electromagnetism, through the n — p mass difference, nuclear 
binding energies, and Coulomb-suppressed and radiative nuclear reactions; and 
the strong interaction, throughout nuclear physics. 

Many studies of the effects on BBN of varying one or more parameters in 
particle physics have been performed in recent years [271 IMl UHl EQl [311 [32l [331 
[34ll35l[36l[37l[38l[39lllQl[4l]; see also [T] and references therein. However, the use 
of BBN faces two major theoretical challenges. First, the degeneracy between 
many different variable parameters, to be compared with the small number of 
observable abundances. Only ^He, ^He, deuterium and ^Li have currently been 
measured to a level of accuracy which gives some prospect of comparison with 
theory, and of these ^He is subject to very large uncertainties in extrapolating 
back to the primordial abundance. We also consider ^Li, since it has been 
suggested that the stellar isotope has been measured [78] and its primordial 
abundance can be predicted, although at such a low level that the primordial 
value may very easily be swamped by ^Li generation through astrophysical 
events or energetic late particle decay. 

Secondly, theoretical uncertainty concerning the way in which the QCD 
parameters, in particular quark masses, affect nuclear forces, and thus nuclear 
binding energies and cross-sections. A derivation of nuclear forces from first 
principles, apart from the long-range attraction attributable to pion exchange, 
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is lackingo Instead, various types of effective tlieories liave been used in order to 
fit measured nuclear properties. These can yield a correct dependence on QCD 
parameters only to the extent that they reflect the underlying physics. There 
are also milder uncertainties in the dependence of certain nuclear reactions on 
a, but in general the role of electromagnetism is well understood. 

In this work we present a unified and systematic approach to these chal- 
lenges, in two respects. First, we consider each variation in particle physics or 
"fundamental" parameters independently. Thus at linear order we can allow for 
any theoretical scenario in which the variations of these parameters are subject 
to some kind of unified relation. Second, we identify which nuclear proper- 
ties and reactions have significant influence on the variation of the primordial 
abundances. This is done by varying every relevant nuclear binding energy and 
cross-section in a reaction integration code, based on the Wagoner /Kawano 
code [l2l 031 SH Us] with updated numerical techniques, and noting the devia- 
tion of output abundances away from the unperturbed values. (Changes to the 
code are discussed in more detail in Appendix A.) This response to the nuclear 
parameters has then to be related to the variation of fundamental parameters 
via nuclear theory. 

The results indicate where the sources of most theoretical uncertainty cur- 
rently are, and suggest where improvements in nuclear theory are most needed. 
Our approach allows us to disentangle these problems. Our study of the de- 
pendence on nuclear parameters can be used as a basis, even if new theoretical 
insights modify the relations between nuclear physics and particle physics pa- 
rameters used in this paper. The connection between the different levels of 
physical understanding arises by simple matrix multiplication of "response ma- 
trices" in a linear theory. 



2 Method 

We consider the set of primordial abundances Ya with a = (D, ^He, ^He, ^Li, 
''Li) and study its dependence on the variation of a set of nuclear physics pa- 
rameters Xi. Here the index i denotes different particle masses such as or 
TTT-e and nuclear binding energies, as well as the neutron lifetime r„ and cou- 
plings and mass scales such as the fine structure constant and Newton constant 
Gn that enter the calculation of nuclear abundances. Our central quantity is 
the response matrix C with matrix elements |3T] 



''"'"ainxr 

It indicates the leading linear dependence for small deviations of the abundances 
about the values obtained given the nuclear parameters inferred from present 
laboratory experiments. The matrix C is extracted by varying the quantities 
Xi independently in the BBN code, a procedure which includes variation of 
the reaction cross-sections and rates that have a physical dependence on Xi. 



■^Recent efforts in lattice QCD [55] have been encouraging. 
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All variations in parameters are taken to be small, such that all necessary 
information can indeed be extracted from the response matrix. 

The second step involves the relation between a set of Standard Model 
parameters Gk and the nuclear physics parameters Xi encoded in a second 
response matrix F with entries 

= (2) 

This second step requires theoretical assumptions and contains, at present, 
substantial uncertainties, which we discuss in Section HI 

The variation of abundances with respect to the fundamental parameters 
Gk is expressed by the "fundamental response matrix" R with elements rak'- 

The matrix R is obtained from C and F by simple matrix multiplication: 

R = CF. (4) 



2.1 Variation of dimensionless and dimensionful parameters 

It is sometimes overlooked that the variation of a dimensionful quantity is not 
physically well-defined: actually, one dimensionful quantity can only be mea- 
sured by comparison with another. The fractional variations of dimensionful 
quantities are thus completely dependent on the choice of units; only dimen- 
sionless ratios are measurable. 

In practice, one may choose a reference frame where some given mass scale is 
kept fixed. Popular frames are the Einstein frame where the Planck mass is kept 
fixed, or the Jordan frame where some particular particle physics scale is kept 
fixed. One may change from one frame to another by a Weyl transformation 
of the metric — the time variation of dimensionless couplings and mass ratios 
is independent of the frame ^ |6jll In this paper we use a frame where the 
QCD invariant scale Ac (sometimes written as Aqcd) is kept constant. This 
is convenient for dealing with nuclear reactions, where the masses and energy 
scales are mainly determined by the strong interactions. Thus the variations 
of dimensionful parameters include implicitly some power of Ac. For example, 
if we take the electron mass mg as a parameter we are implicitly considering 
a variation of mg/Ac. The implications of varying dimensionful parameters in 
QCD and nuclear physics are discussed further in Appendix lA. 31 



3 Nuclear response matrix 

We first establish the matrix C by a systematic variation of parameters in 
the BBN code. For this purpose we have modified the code as described in 
Appendix A, which also contains the details of our treatment of reaction rates. 

^This supposes that the time interval is appropriately rescaled. However this is not an 
issue for BBN, since we do not consider the time derivative but rather the absolute variation 
between BBN and the present time. 
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3.1 Nuclear parameters and response matrix 

We vary with respect to the following thirteen "nuclear" parameters Xi in the 
BBN code: 

• Gravitational constant Gm 

• Neutron lifetime t„ 

• Fine structure constant a 

• Electron mass mg 

• Average nucleon mass mjv = (m„ + nip)/ 2 

• Neutron-proton mass difference Qn = mn — frip 

• Binding energies of D, T, ^He, ^He, ^Li, ''Li, '^Be 

Clearly to vary each binding energy independently is unphysical, however our 
purpose is to determine the leading (linear) variation of all abundances, with 
respect to the code's input parameters, once and for all. Then given any specific 
theoretical model, we can construct a linear combination of variations of Xi to 
account for any variation of a fundamental parameter. 

We should also consider the possibility that some states could change be- 
tween bound to unbound, or stable and unstable, under a variation of some 
fundamental parameters. As discussed in Section [321 we do not consider vary- 
ing the binding energies Bi oi A <7 nuclei past the point where the Q- value of 
any reaction relevant to standard BBN changes sign. If some Q-values approach 
zero we already find large variations in abundances. However, we should also 
check whether some reactions that are unimportant in standard BBN become 
relevant under a small change in parameters. If this were so then the linear 
approximation would quickly become inaccurate. 

One could expect substantial changes in the final abundances of BBN if 
either the dineutron were bound, or the *Be nucleus were stable, at the time 
of nucleosynthesis. Such effects would not be seen in a purely linear analysis 
expanding about the present-day values. We discuss the dineutron and ^Be 
further in Section 14.41 We find that the dineutron binding cannot produce 
significant effects, given the other observational bounds on the variation of 
couplings. The requirement that '^Be should remain unstable and short-lived 
during BBN gives a one-sided bound on the variation of a, but only if light 
quark masses are held constant relative to Ac. 

Our results for the nuclear response matrix are shown in Table [TJ The first 
thirteen rows constitute the transposed nuclear response matrix C'^ . We also 
quote the dependence of the abundances on rj in the last row. Values are quoted 
to 2 d. p. or to 2 sig. fig. when the magnitude exceeds 1. Below we give a few 
comments concerning specific parameters Xj. 
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Table 1: Response matrix C, dependence of abundances on nuclear parameters. 

Gravitational constant: A variation in Gn or equivalently the Planck mass 
Mp, relative to Ac, affects the Hubble expansion. It is thus equivalent to chang- 
ing the energy density of the Universe by a constant factor. During BBN this 
is also equivalent to adding or subtracting a number of relativistic species with 
radiation-like equation of state, the so-called "effective number of neutrinos"; 
and also to the presence of a 'tracker' form of early dark energy [4]. 

Neutron lifetime / Gp: The neutron lifetime appears in normalising the 
n ^ p reaction rates and the free neutron decay. This is essentially the only 
important weak interaction, hence a variation in r„ is equivalent to a variation 
of GfE More precisely, for fixed mjy, me and Qn and Qa/dv (the ratio of 
nucleon axial vector / vector couplings), there is a one-to-one relation between 
r„ and Gp such that we can choose to work with one or the other parameter 
equivalently, using the mapping Alnr„ = —21S.\\iGf- Our results for variation 
of Tn are then consistent with [271 ED [33] . 

Fine structure constant: Charged particle reactions incorporate an explicit 
a-dependence in the Gamow factor, reflecting the leading effects of tunneling 
through the Coulomb barrier. Coulomb corrections also affect the n <-> p weak 
reaction rates. As noted in [36], radiative reactions including npd'y clearly have 
an additional dependence on a. Reactions with charged particles in the final 
state also experience Coulomb suppression, though generally much less signifi- 
cant since the outgoing momenta are larger. For Q > E the final state Coulomb 
suppression is negligible; the only cases where it may be worth considering are 
^He(n,p)3H and ^Be(n,p)^Li, with Q-values of 0.76 MeV and 1.64 MeV respec- 

*The variation of t„ is also important in accounting for uncertainties in the laboratory 
value. 
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tively. Such effects were estimated by Nollett and Lopez using the Coulomb 
wavefunctions for orbital angular momenta I = 0,1,2 and appear to be small, 
the strongest dependence being (av) oc a~^'^. 

Our calculation of the a-dependence as a nuclear parameter does not in- 
clude the effect of varying a on nucleon masses and nuclear binding energies; 
these are varied independently. The a-dependence of binding energies will be 
incorporated at the second step through the response matrix F. On this second 
level we treat a as a fundamental coupling and incorporate the a-dependence 
of the nuclear physics parameter Xi. 

Neutron, proton and electron masses: We explicitly calculate the de- 
pendence of the n ^ p weak reactions on the electron mass nie, and divide up 
the nucleon masses into an isospin-invariant part raj^ = (nip + m„)/2 and an 
isospin- violating part Qjy = nip — nin ■ The dependence of final abundances on 
these nuclear parameters appears somewhat counterintuitive, since at this level 
we vary them with the condition that r„ should remain constant. Thus effec- 
tively the variations in nie and Qn are accompanied by compensating variations 
in the Fermi constant. When this is taken into account, our results for variation 
of Qn are consistent with those of \27\ \31 \ (However, our treatment of rUe 
as a nuclear parameter differs from both [31] and ^27j hence the results are not 
directly comparable.) 

The variation with respect to m^v = (m-n + mp)/2 is mainly due to the cal- 
culation of thermal reaction rates {av) at given temperature, since the relation 
between energy and velocity is affected; thus abundances which depend strongly 
on reaction rates will feel this variation. Again as with binding energies, we 
cannot at present tell to what extent reaction matrix elements depend on mjv 
(with the exception of npdj). However, for any choice of variations in funda- 
mental parameters, we will see that the variation oim]\j is very small compared 
to other nuclear parameters, since we have fixed Ac to be constant. 

3.2 Variation of binding energies 

Table [J indicates that by far the largest sensitivity of abundances to nuclear 
parameters involves the variation of ^He, ^Li and ^Be binding energies. This 
can easily be understood since the reaction rate of ^He(a,7)^Be with Q-value 
1.59 MeV is very sensitive to changes in these (numerically large) binding en- 
ergies. In order to have a possibly more intuitive idea of the sensitivity of Ya 
to binding energies we display in Table [2] the dependence of abundances on 
the binding energy per nucleon hi = Bi/Ai, that is, dlnY/dbi, since it is more 
meaningful to compare values of bi between nuclei of different A. We see that 
the dependence of the lithium abundances on the ^He, ^Li and ''Be binding 
energies is still dominant. 

Each reaction Q-value is determined by the masses of reactants and prod- 
ucts: the relevant parameters are the binding energies of nuclei up to A = 7. 
The Q-value of each reaction affects the abundances via the reverse thermal 
reaction rate relative to the forward rate, and via phase space and radiative 
emission factors in the reaction cross-sections. 
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Table 2: Dependence of abundances on binding energy per nucleon, in units of 
MeV-^ 



The reverse reaction rate is simply related to the forward rate via statistical 
factors, due to time reversal invariance (see for example [S]): the relevant 
dependence is 

^"oce-Q/^. (5) 



12^34 

The Q-dependence of radiative capture reactions (assuming a dominant electric 
dipole) is 

a{E) (xE?^r^{Q + Ef (6) 

whereas for 2^2 inelastic scattering or transfer reactions the dependence is 

a{E) oc f3 ^ {Q + E)^/^ (7) 

where /? is the outgoing channel velocity. In the current treatment we assume 
-C Q at relevant temperatures, and simply scale rates by the appropriate 
power of Q. Clearly this breaks down when Q approaches zero, and we have 
not considered varying any binding energy to the point where this happens. A 
more accurate treatment would involve applying the phase space dependences 
directly to the cross-sections, for example in the S-factor description of charged 
particle reactions, which involves an expansion in E; the dependence on (Q + E) 
can then be applied order by order. 

For resonances, which appear in reaction rates as terms varying with e"^''/^, 
we scale their contributions by the appropriate power of iQ+Er)n The question 
of whether E,. should be varied is still open. The minimal assumption is that 
the mass-energy of the resonance varies with the mass-energy in the incoming 
channel, thus Er would be constant to first approximation. 

A variation in binding energies can have two kinds of effect. The Q value 
can change the time when a reaction drops out of equilibrium, for instance the 
n+p (i+7 reaction. Or it can change the absolute rate of a reaction, and thus 
the production rate of a given species, for example the ^Be-producing reaction 
whose cross-section varies with Q^. 

Whether the reaction matrix elements have a dependence on the binding 
energies and on Q is in general not clear because there is no systematic effec- 
tive theory for multi-nucleon reactions. The exception is the npdj reaction. 



^The resulting variations of abundances are indistinguishable from the result of scaling by 
a power of Q, within our uncertainties. 
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for which we have implemented a nuclear effective theory result [51j where de- 
pendence on Bj) is explicit. Other parameters of this effective theory may also 
have some quark mass dependence, and there are ongoing efforts to connect it 
to chiral perturbation theory and thus to QCD parameters, specifically mq/Ac 



It has been suggested (see [S^) that the cross-sections of other reactions in- 
volving deuterium, specifically when D appears in the ingoing channel, will have 
a significant additional dependence on Bj). To first order, the scattering length 

— 1/2 

varies as Bj^ , and cross-sections in the low-energy limit may be expected 
to vary with B^^, in addition to previously discussed effects. In Tabled] we 
present values for the dependence on Bd, assuming this additional dependence 
not to be present. When, however, a o" oc 1/Bd dependence is implemented for 
the 7)^He, d{d,n)^I{e and d{d,p)t reactions, we obtain 

f/"^" (D, ^He, ^He, ^Li, ^Li) = (-1.5, -2.4, 0.66, -5.6, 7.5). (8) 

The resulting uncertainties in reaction rates have a subleading effect on abun- 
dances, compared to other effects of varying B^. To be consistent we should 
similarly consider the effects of all other binding energies on reaction cross- 
sections beyond the kinematic factor of Eqs. (j6l7p : however in any case it is 
not clear whether the scattering length is the correct parameter to consider. 
Therefore we will take the case without this a oc l/Bj) dependence as our fi- 
nal result; the other case serves as an illustration of possible further effects of 
binding energies on reaction rates. 

Our results for the dependence on Bd are consistent with those of [29^ 
[3H (allowing for the different treatments of D-destroying reactions) but 
differ from [27| and [32] for the dependence of the D and ^Li abundances. The 
discrepancy with respect to those works arises from a different treatment of the 
npd'j reaction. It is also not clear to us if the effect of changing Bd on other 
reaction Q-values and rates was included. 



3.3 Nuclear rates 

In order to estimate which reactions are more or less important in the variation 
of the final abundances, we varied each thermal averaged cross-section {av) 
by a temperature-independent factor, preserving the relation between forward 
and reverse rates. The aim is to diagnose which reactions one should focus on 
in discussing the sensitivity of observed abundances to variations of couplings. 
The n <^ p weak interactions influence every abundance nontrivially (we will 
treat them analytically); in addition, based on our results, we designated cer- 
tain other reactions as important. These are given in Table [sl^ We judge a 
reaction cross-section to be important if the dependence of any abundance on 
any given reaction cross-section dlnYa/d\n{av)i is more than 0.1. In every case 
the dependences are order unity or smaller, and for many "important" reactions 

®As usual in BBN simulations, the slow /3-decays of tritium and ^Be are accounted for by 
adding on the T and ^Be abundances to ^He and ^Li respectively at the end of the run, when 
other nuclear reactions have frozen out. 
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Table 3: Leading dependence of abundances on thermal averaged cross-sections 
d\n.Ya/ d\n.{av)i for important reactions (1 d. p.) 

only a few abundances are significantly affected. A reaction cross-section can be 
unimportant (for our purposes) for one of two reasons: either it is so small that 
the reaction is irrelevant, or the reaction is so rapid, compared to the Hubble 
rate and to other slower reactions, that the rates of change of abundances are 
almost independent of the cross-section. Hence our list of "important" reactions 
differs from [48j . 

Note also that the ^He abundance does not depend on any nuclear reaction 
cross-section (apart from n ^ p). Furthermore, the reactions involving ^Li are 
not important for any other measurable primordial abundance. 

In implementing the variations of nuclear parameters as displayed in Ta- 
bled! we do not directly use the effects of a temperature-independent variation 
of integrated cross-section {av) given in Table El Since variations of Xi in gen- 
eral result in temperature-dependent variations of reaction rates, we implement 
them directly within the integration code. 

4 Relations to fundamental parameters 

The next task is to connect these nuclear parameters to fundamental param- 
eters Gfc at a higher energy scale. We consider the following six fundamental 
parameters G^: 

• Gravitational constant Gm 

• Fine structure constant a 

• Electron mass mg 

• Light quark mass difference 5q = rrid — m„ 

• Averaged light quark mass rh = [rrid + m„)/2 oc 

• Higgs v.e.v. {(p). 
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Table 4: Response matrix F, dependence of nuclear parameters Xi on funda- 
mental parameters Gj. 

An additional parameter is the strange quark mass m^. We have omitted it 
from our list because the present theoretical uncertainties of how it influences 
the nuclear parameters are too high. 

The dependence of the fundamental parameters G^ on the nuclear param- 
eters Xi is encoded in the matrix F defined in Eq. ([2]): our estimates of F are 
shown in Table HI We now discuss how to derive the entries of this matrix. 

The gravitational constant enters as before, corresponding to the (1, 1) entry 
of unity in Table HI The fine structure constant influences the abundances in 
two ways. First we have the direct influence of its variation while keeping 
other nuclear parameters fixed. This is accounted for by the value 1 in the 
(2, 2) element. Secondly, it influences the nucleon masses and nuclear binding 
energies. This yields the other elements in the second column of Table HI As a 
nuclear parameter a does not depend on (0), me, etc., as reflected in the second 
row of Table HI 

The a-dependence of the neutron lifetime (91nT„ = 3.86 5 In a enters via 
the n-p mass difference. The full dependence of on fundamental parameters 
can be found in [31]. Similarly as for a, the electron mass as a fundamental 
parameter has an effect on the neutron lifetime, dluTn = 1.52 9 In mg, and thus 
also on n p reaction rates. 

Elementary u and d quark masses appear as fundamental parameters. They 
may be divided up as {mu+md)/2, which influences the pion mass, which in turn 
determines the behaviour of nuclear forces; and — m„, which influences the 
neutron-proton mass difference. The average light quark mass also affects the 
nucleon mass tjin via the so-called sigma term; this holds also for the strange 
quark mass through the "strangeness content" of the nucleon. The effect of 
varying pion (i.e. quark) mass on nuclear binding energies and reaction cross- 
sections is in general not known. However, for the binding energy of deuterium 
and the npd'j reaction the dependence has been found to some approximation 
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using effective tlieories. 

At tliis stage we vary the Higgs v.e.v. independently of the elementary 
fermion masses, thus it only influences weak reactions via Gp: the relevant 
reactions are neutron decay and p. 



4.1 Fine structure constant 

The dependence on a of the masses of (composite) particles and nuclear binding 
energies is found by estimates of electromagnetic self-energy or binding energy, 
for example [58] for the proton-neutron mass difference, and the semi-empirical 
mass formula for nuclei. Note that the fractional variation of ttt-tv with a is 
negligibly small. More precise estimates for nuclear binding energies have been 
made using realistic models of nuclear forces (see [59] and [60]) and similar 
values appear in [36]. Variation in a as a fundamental parameter leads to 
variations of nuclear parameters: 

A ln(r„, Qtv, -Bd, -Bt, -Bshg, B^^ie, -BeLi, ^7Li, ^7Bc) = 
(3.86, -0.59, -0.0081, -0.0047, -0.093, -0.030, -0.054, -0.046, -0.088) A In a. 

(9) 

Then the resulting variations of abundances are 

A ln(yD, l3He, l4He, IfiLi, l7Li) = (3.6, 0.95, 1.9, 6.6, -11)A In a (10) 

These results are similar to those of |40j and ^36j. The "^Li abundance here 
is more sensitive to a than the Nollett &: Lopez estimate [36]. This appears 
due to their use of a cluster model for the ^He(a,7)'^Be reaction, whereas our 
treatment of this reaction simply uses the Gamow factor. 



4.2 Electron mass 

The electron mass affects the neutron lifetime via 

d In T„, 



9 In 



1.52. 



Adding this variation to the effects already calculated gives, for the variation 
of me as fundamental parameter: 

Aln(yD,^3He,nHe,l6Li,lVLi) = (0.46, 0.21, 0.40, 0.97, -0.17) A In r^e. (11) 

Our result for the variation of me as fundamental parameter is very similar to 
the semi-analytic result of |31j . 



4.3 Pion and light quark masses 

The effect of light quark masses on the nucleon masses can be found by con- 
sidering low-energy hadron physics (see for example [581 [6T] ) . We consider an 
average light quark mass m = (m^j + mrf)/2, the mass difference 5q = m^ — m^ 
and the strange mass m^ . Note that the m^ dependence of hadronic physics is 
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still quite unclear, for example the strangeness content of the nucleon is subject 
to at least 50% uncertainty. The nucleon mass gets a contribution from nonzero 
quark masses 

rriN = m^^(Ac) + fT,rAf(m.) (12) 
where is the mass in the chiral limit and 

TTl — 

(^ttN = 7, {p\uu + dd\p). (13) 

For {p\uu + dd\p) depending only on Ac, one has din aj^jsf /dlnrh = 1 and 
therefore 

= 0.048 (14) 

a In m m at 

and for the strange quark 

^^2(Mss}p)^ ^ §^ = !I^^^ 0.12 ±0.12 (15) 
{p\uu + dd\p) omms m IrUn 

given nis/rh ~ 25 and y = 0.2 ± 0.2 |56j . 

The strange quark mass is close to the nonperturbative scale Ac. For this 
reason a systematic or physically meaningful treatment of the dependence of 
nuclear quantities on has not been possible so far. Thus our results have 
the caveat that variation in m^/Ac is not yet accounted for. 

The dependence of the nucleon mass difference Q on quark masses was 
estimated in |58j : in units where Ac is constant, we have 

AQiv ^ (-0.76Alna + 2.05Aln(5„)MeV ^ t!^? - 1.59. (16) 

A m bq 

Recent lattice QCD studies with dynamical quarks [62] have calculated the 
dependence of the nucleon mass splitting on 5q: the result is consistent with 
the estimate we adopt. 

The pion mass is crucial for nuclear forces and its leading dependence on rh 
follows from chiral perturbation theory as 

= m{uu + dd)f~'^. (17) 

Here we assume the leading order where and {uu + dd) depend only on Ac, 
therefore AlnmTr — ^Alnm. 

Static properties of nuclei, and importantly for BBN, nuclear binding ener- 
gies, depend strongly on the pion mass, which determines the range of attractive 
nuclear forces. Quantum Monte Carlo calculations have been performed with 
realistic nuclear potentials and accurately reproduce many experimental prop- 
erties [60|. i59j . One-pion exchange and two-pion exchange are dominant contri- 
butions within the expectation values of the two- and three-nucleon potentials 
respectively. Currently such studies have not been extended to determine the 
functional dependence of binding energies on the pion mass in general. This de- 
pendence would in any case have uncertainties due to subleading effects of pion 
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mass (or equivalently light quark masses) on other terms in the nucleon-nucleon 
potential [65] , 

However, the dependence of the deuteron binding energy on the pion mass 
has been extensively studied within low-energy effective theory \64:\ I65j : the 
result may be expressed as 

r 

AIuBd = rAlnrriT^ = —Alnm (18) 

for small variations about the current value [35], with — 10 < r < —60 We 
will also take this dependence as a guide for the likely pion mass dependence 
of other binding energies. Although the size of the deuteron binding appears 
due to an accidental cancellation between attractive and repulsive forces, its 
derivative with respect to (which is just Bd/ttIt^ times r) is not expected 
to be subject to any cancellation. We also expect that the pion contribution 
to the total binding energy should increase with the number of nucleons; a 
proportionality to {A — 1) seems reasonable to obtain correct scaling at both 
small and large A. Hence to estimate the effect of pion mass on the binding 
energy of a nucleus Bi we set 

1^ = MA, - l)^r -0.m{A^ - 1) (19) 

taking r ~ —8. The numerical constants fi are expected to be of order unity, 
but will differ between light nuclei due to peculiarities of the shell structure, etc. 
Our normalization corresponds to /d = 1. Then the nontrivial dependences of 
nuclear parameters on m are 

A ln{B-[), B^, -BsHe, B^e, Bqu, Bru, B^BcrnN) — 
(0.5r, 0.26/Tr, 0.29/3Her, 0.12/4Her, O.Uhur, 0.17 fmr, 0.18/7Bcr, 0.048) A In m. 

(20) 

For the rh dependence of abundances due to the variation of binding energies 
we then have 



d In Ya 



dlnm 



r^^ iA-l)Bj,dlnY, 



Taking account also of the small effect of rh on the nucleon mass mjy, the 
resulting dependence of abundances on m is 



d\nm. 
dlnm 

dlnm 
dlnVeu 

dlnm 
dlnYju 

dlnm 



11 + 0.5/t + 5/3He 
8 + 3/t - 7/3He 

-2.7 

27 + 0.4/t + 7/3He + 55/4He " 96/6Li 

-36 + 5/t + 22/3He + 54/4He + 9/7Li " 115/7Be. (22) 



^Our definition of r differs by a sign from [35] . 
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Even if we consider that some contributions could cancel against one another 
due to the values of the fi, the magnitude of these variations is striking, partic- 
ularly concerning the lithium abundances. To get an idea of the possible effect 
of cancellations, we may set all fi to unity and find the dependences 

Aln(yD,l3He,nHe,l6Li,l7Li) ^ (17, 5, -2.7, -6, -61) A In m. (23) 

One may also consider to what extent varying rh or the pion mass may affect 
reaction cross-sections beyond the npdj reaction. It seems very likely that 
matrix elements would acquire nontrivial dependence on m^; however, since 
the dependence of abundances on reaction cross-sections is relatively mild (see 
Table ^ , the dependence via reaction matrix elements is unlikely to compete 
with the very large effects arising through the variation of binding energies. 

4.4 Stability of dineutron and ®Be 

The effects of a bound dineutron on BBN were studied in ^46j, where it was 
found that the final abundances were essentially unaffected as long as the dineu- 
tron binding energy remained smaller than i?D . An effective field theory analysis 
may also be applied to the binding of the dineutron, which lies in a different 
channel from the deuteron [65]. The sensitivity of the dineutron system to 
the pion mass, or equivalently to m, is found to be comparable to that of the 
deuteron. It is very unlikely that the binding of the dineutron could become 
significant, since the fractional variations of Bd for the range of variations of 
couplings considered in this paper are at most of order 5%, which amounts to 
Ai?D — 0.1 MeV (see Section l5.2p . In earlier work a simple potential model 
was used |39j to investigate the deuteron and dineutron binding, with a similar 
result: the variation in the quark mass parameter rh/Ac required to bind the 
dineutron is much larger than that required to cause even a 100% variation in 
Bd. 

The decay of '^Be to two ^He nuclei has a Q-value of only 0.092 MeV, thus 
it is conceivable that even a very small variation of parameters could result 
in a stable (or long-lived) ^Be at the time of BBN. This would have dramatic 
consequences for the abundances of heavier nuclei since it would then be possible 
to synthesize carbon directly by two-body reactions. Hence we can immediately 
rule out any variation that causes this Q- value to change sign. 

Electromagnetic contributions to the *Be binding are estimated from the 
results of quantum Monte Carlo calculations j£Oj as before. We find A In -BgBe — 
— 0.0588A In a, thus the Q-value varies as 

AQ8^2q ^ (1-60 MeV)A In a. (24) 

Large negative values of Aa are excluded: we find a firm limit Alna > —5.7%, 
under the condition that other fundamental parameters that may affect the 
sign of Q, i.e. quark masses, do not vary with respect to Ac. However, we will 
obtain stronger limits on the variation of a by considering the response of the 
observed BBN abundances: see Section [5.2[ 
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d In Yn Id In Gt 


D 


^He 


''He 


«Li 






0.94 


0.33 


0.36 


1.4 


-0.72 


a 


3.6 


0.95 


1.9 


6.6 


-11 


{^) 


1.6 


0.60 


2.9 


5.5 


1.7 


me 


0.46 


0.21 


0.40 


0.97 


-0.17 




-2.9 


-1.1 


-5.1 


-9.7 


-2.9 


m 


17 


5.0 


-2.7 


-6 


-61 


rj 


-1.6 


-0.57 


0.04 


-1.5 


2.1 



Table 5: Response matrix R, dependence of abundances Yi on fundamental 
parameters Gk 

The Q-value for ^Be decay results from a nearly exact cancellation between 
®Be and ^He binding energies, thus its dependence on light quark masses is sub- 
ject to very large theoretical uncertainty. We again estimate the m-dependence 
in terms of the deuteron binding and find AlnSgBe = — l-l/gBo Alum. Thus 
for the Q- value, 

AQ ~ ([9 + 62(/8Bc - 1) - 53(/4He - 1)] MeV)Alnm. (25) 

The expression inside brackets depends very strongly on the unknown /j factors: 
its likely magnitude is 10-50 MeV barring cancellations. We would then have 
a one-sided bound on variation of rh at the 1% level or better from stability of 
®Be. At present we do not know the sign of the prefactor and the use of ^Be to 
place useful bounds on the variation of QCD parameters must be left to future 
work. 



5 Dependence of abundances on fundamental pa- 
rameters 

5.1 Fundamental response matrix 

We next combine the nuclear response matrix, Table [U with the relations be- 
tween nuclear and fundamental parameters (Table H]), according to Eq. (jj]). 
Table [5] shows the resulting dependences of abundances on fundamental param- 
eters, as encoded in the matrix R. This table is our central result. 

In treating the m-dependences, which arise from the nuclear binding energies 
with their uncertain values of /«, we have given the values which arise when 
setting all fi to unity. Alternatively, if all the are of order unity but the term 
with the largest prefactor dominates, we would obtairH 

|j^(D,3He,4He,6Li,7Li) = (11, 8.4, -2.7, -96f6Li, -115f7Be). (26) 
a mm 

The dependence on Gn is consistent with the results of [331 EI] and p9| 
I30j . once one translates from units where Gn is constant to ours where Ac is 
constant. 

*If, on the contrary, there is substantial canceUation then the dependences of the deuterium, 
^He and lithium abundances on m may be smaller. 
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-19% 


< 


AlnGjv 


< 


+10% 


-3.6% 


< 


A In a 


< 


+1.9% 


-2.3% 


< 


A ln(0) 


< 


+1.2% 


-17% 


< 


A InTTI-e 


< 


+9.0% 


-0.7% 


< 


Aln6q 


< 


+1.3% 


-1.3% 


< 


A Inm 


< 


+1.7% 



Table 6: Allowed individual variations {2a or "conservative allowable range") 
of fundamental couplings 



The ^He dependence was previously calculated in [3T] by semi-analytic 
methods: our results for the dependence on fundamental parameters are similar. 

The {(p) dependence may be compared to the result of [30] for the variation 
of Gp oc (using a semi-analytic method): the results for ^He, ^Li and ^Li 

are similar, but much smaller variations of D and ^He are obtained in [30j . 



5.2 Bounds on separate variations of fundamental couplings 

The first application of the result is in setting bounds on the variation of each 
fundamental parameter considered separately, under the assumption that only 
one parameter varies at once. We may consider three observational determi- 
nations of primordial abundances (see Appendix [B|): deuterium, ^He and ^Li. 
However, the observed ''Li abundance deviates by a factor two to three from the 
value predicted by standard BBN theory (SBBN), and systematic uncertainties 
related to stellar evolution exist [79]. Thus, we use the former two, D and ^He, 
to constrain the allowed variations of the fundamental constants individually. 
For deuterium we take 2a limits; for ^He we consider instead the "conservative 
allowable range" of [69] . The resulting constraints are given in Table [H 



5.3 Sensitivity matrix 

The relative precision of the observational determination of primordial abun- 
dances is best for ^He and somewhat poorer for D and ''Li. This can be taken 
into account by defining a "sensitivity matrix" S with elements Sak- 

-I dYa ia, th /„f7\ 
^ak = <7a a^ - ^"^.th (27) 

where the subscript 'th' denotes a quantity evaluated about the values obtained 
in SBBN, and the approximation follows from taking a linear dependence. Here 
aa is the \a error of the observational determination of the primordial abun- 
dance Ya- thus Ya/aa is a measure of the precision with which a given abun- 
dance is known. The current situation of theory and observation is summarized 
in Appendix [BJ We adopt SBBN theoretical values 

yD,th = 2.61x10-5 
nHe,th = 0.2478 

l7Li,th = 4.5x10-1° (28) 
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a-'dYjd\nGk 


D 




'Li 


Gn 


6.1 


9.9 


-6.5 


a 


24 


52 


-100 


{^) 


10 


79 


15 




3.0 


11 


-1.5 




-19 


-140 


-26 


m 


110 


-74 


-550 



Table 7: Sensitivity matrix S 



and la observational errors 

fJD = 0.4 X 10"^ 

0-4He = 0.009 

a^u = 0.5 X 10"^°. (29) 

We then obtain the sensitivity matrix S given in table [71 

The matrix elements Sak have a simple interpretation. If we perform a 
variation of a given fundamental coupling Gk by 1% while keeping all other 
fundamental couplings fixed, the variation of Ya/aa is given by Sa^/lOO. The 
range of variation of each Gk corresponding to the la (etc.) range of any ob- 
served abundance can easily be read off. Thus for example the currently allowed 
la range of Yme [69] corresponds to a variation in a of 1/52 ~ 1.9% each side of 
a central value, or 3.8% in total. Considering individual variation of couplings 
we find that at present the '^He abundance is the most sensitive measurement 
of the couplings Gn, 0, "i-e and 6^; on the other hand a and the average light 
quark mass m are most limited by ''Li. 

5.4 Inverse sensitivity matrix 

A quantity like the inverse of the sensitivity matrix would also be useful. With 

6lnGk = tka<T-^6Ya, (30) 

we could immediately infer the fractional change in fundamental couplings cor- 
responding to a deviation of one abundance Ya away from the standard pre- 
dicted value, while keeping the other abundances fixed. For example, a 3a 
change in ^Li, i.e. 6YYhi/o'7u = 3, would be correlated with a simultaneous frac- 
tional change 3tka in the couplings Gk- Of course, if we have more couplings 
than observed abundances, the matrix T with elements t^a is not uniquely 
defined. We may, however, restrict the number of independently varying fun- 
damental couplings, either by keeping some couplings fixed or by assuming 
relations between the Gk, for example motivated by grand unification. If S is 
reduced in this way to a n x n matrix one simply has 

T = S~\ (31) 
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a-'dYjd\nGk 


D 


^He 


''Li 


Gn 


6.1 


9.9 


-6.5 


a 


24 


52 


-100 


(0) 


100 


-125 


-560 





D 


^He 




Gn 


0.24 


-0.015 


0.048 


a 


-0.036 


0.016 


-0.010 


{^) 


0.0038 


-0.0027 


-0.00048 



Table 8: Sensitivity matrix S and its inverse T, given constant Yukawa couplings 



As a simple example, we may discard the variation of the Yukawa couplings 
such that 

A In rfi = A In 5g = A In m-g = A ln((/)) . 

For this case we show the matrices S (transposed) and T in Table El A decrease 
of 6(T for ^Li brings its predicted abundance from 4.5 x 10^^'^ to 1.5 x 10~^'^. 
With other abundances fixed, this could (in linear approximation) be realized 
by the variations 

Aln(GN,a, ((/>)) = (-29%, +6.1%, +0.29%). (32) 

The required changes in Gn and a are relatively large: our linear treatment may 
still be valid here, but in general this should be checked, which we do explicitly 
for specific unified models in Section 16. 5i Alternatively, for a f>a decrease in 
''Li from the standard predicted value and a la increase in deuterium, bringing 
the predicted Yd to 3.0 x 10~^, we would require variations of fundamental 
parameters 

A ln(GN, a, (</>)) = (0, +1.4%, +0.62%), (33) 

clearly well within the linear regime with respect to variation of fundamental 
parameters. 



6 Unified models 

It is of interest to consider unified scenarios where the variations of fundamen- 
tal couplings satisfy relations that reduce the number of free parameters. In 
the simplest case every variation of a parameter G^ is determined by a single 
underlying degree of freedom. The variations can then be written as a vector: 

AlnGfc = 4A(^ (34) 

where (/J is a (dimensionless) field which gives rise to the variation and are a 
set of numbers characterising a particular unified model. We then obtain 

Alny, = (GF),fc4A(^, (35) 

where we may also eliminate A99 in favour of the variation of some observable 
parameter. 

We will consider four simple possibilities. First, that the strength of gravita- 
tion varies, but all other scales and couplings of particle physics are unchanged. 
This corresponds to a violation of the strong equivalence principle, while the 
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weak equivalence principle is preserved. Thus has a single nonzero entry 
corresponding to Gn and we find 

Aln(yD,l3He,nHe,l6Li,nLi) ^ (0.94, 0.33, 0.36, 1.4, -0.72) A IuCGnA^) (36) 

restoring the implicit dependence on the QCD scale. 

In the remaining examples, we consider a grand unified theory with unified 
coupling ax, broken at the scale Mx to the Standard Model symmetry. The 
observable couplings of QCD and electromagnetism are then related to ax via 
renormalization group (RG) flow and electroweak symmetry breaking. Thus one 
can find a relation between the variation of a and that of h.c/Mx, depending 
on the gauge group and matter content of the theory [66l [39l |T0] . 

We also need to specify the behaviour of the ratios of energy scales Mx /Mp, 

— 1/2 

{(j))/Mp, me,g/Mp, where Mp is the Planck mass proportional to . In all 

unified scenarios we will take the Planck mass fixed relative to the unification 
scale, thus A{Mp/Mx) = 0. For simplicity we take the Yukawa couplings to be 
constant, thus the electron and quark masses are proportional to ((/>). We also 
define an exponent 7 which relates the variation of {(J)) with respect to Mx to 
the variation of Ac/Mx as 

Mx \Mx) ^ ' 

6.1 Fixed ratio of weak and strong scales 

In the first unified scenario "GUTl" we take all low-energy mass scales of 
particle physics to be proportional to Ac, thus 7 = 1. Then using the non- 
supersymmetric GUT relations discussed in [311 IHl [10] we have 

22a 

A In a = Alnax ^ 0.92Alnax, (38) 

7ax 

Aln— = — Alna ~ 39Alna (39) 
Mx 11a 

taking ax — 1/40. These relations include the effects of varying (relative 
to Mx) charged particle masses, or "thresholds", on the variation of gauge 
couplings H Then we have 

Aln(GN,a, (0),me,5g,m) ~ (78, 1, 0, 0, 0, 0)A In a. (40) 

In this case the variations of abundances are not subject to the theoretical 
uncertainty of varying m^/Ac. We obtain 

Aln(yD,l3He,nHe,l6Li,l7Li) ^ (77, 27, 30, 120, -68) A In a (41) 

where for definiteness we have taken all /j to unity in Eq. ()19p . Note that when 
the variations are reexpressed in terms of A In Gn as 

Aln(yD,1^3He,l4He,l6Li,l7Li) ^ (0.99, 0.34, 0.38, 1.5, -0.87) A In Gn (42) 



In supersymmetric unified theories the efi'ect of varying thresholds (charged particle 
masses) on the gauge couplings may be much larger |67) . 
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the result is very similar to the first scenario where only Gn is varying. The 
scenarios only differ by a variation Aa/a = (1/78)AGn/Gn, hence we do not 
plot separately the first scenario of varying only Gn • 



6.2 Fixed weak scale and varying strong scale 

In the second unified scenario "GUT2" we consider that the Higgs v.e.v. and 
elementary fermion masses are all proportional to the unification scale, thus 
A(Mx/Mp, ((/))/Mp,me,g/Mp) = 0, or equivalently 7 = 0. Then, after con- 
verting to QCD units, the mass scales Mp, {(f)) and rrie^q will vary inversely to 
K/Mx- We find [3T] 

8a 

Alna(Mvi/) = A In ax, (43) 

Due to the effect of the three light quarks whose effect on the running of a(/u) 
is cut off at /i ~ Ac, the variation of the fine structure constant is as follows 

-Alna = — Alna(MvK) ( 1 + T^> ^Q- 1 (45) 



w 



where i runs over three colors of u, d, and s quark, thus = 2. Then the 

variations of fundamental couplings are related as 

Aln-^ = Aln^ = -^Alna ~ 32.3Alna (46) 

Mx {(t>) 40a ^ ^ 

and we have 

Aln(GN,a, ((/>),me,5g,m) ~ (64.5, 1, -32.3, -32.3, -32.3, -32.3) A In a. (47) 

We then obtain variations of abundances 

Aln(yD,>3He,nHe,^6Li,^7Li) ^ (-450, -130, 170, 380, 1960) A In a (48) 

Note that this model is subject to additional uncertainty due to the variation 
in the strange quark mass relative to Ac; however it seems unlikely that this 
variation would produce significant cancellations. 



6.3 Varying the weak scale faster than the strong scale 

In the third unified scenario "GUT3" we consider the case when the Higgs v.e.v. 
and fermion masses vary more rapidly (with respect to the unification scale) 
than the QCD scale Ac does: thus 7 > 1. We take 7 = 1.5 and find that the 
variations of fundamental couplings are related as 

Aln(GN,a, ((/>),mc,5g,m) ~ (87, 1, 21.5, 21.5, 21.5, 21.5)A In a. (49) 

The variations of abundances are then 

Aln(yD,^3He,nHe,l6Li,lVLi) ^ (430, 130, -65, -60, -1420) A hi a. (50) 
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6.4 Results 



In Fig. [T] we show the abundance variations given by the three GUT models, as 
a function of the variation of a. Note that we plot only the linear dependence 
of abundances on a, therefore if Alnl^ becomes larger than 1 (as in the case 
of ^Li) the results may be affected by higher order terms. We also show the la 
observational bounds as highlighted regions. Also included in the plot is the 
effect on the standard BBN predictions of varying the baryon-to-photon ratio 
r/ over the 2a range allowed by WMAP 3 year data, 5.7 < 10^°r/ < 6.5. 

It can be seen that in the "GUT2" scenario a reduction of a by about 

0. 025% (i.e. a fractional variation of —2.5 x 10"'*) would bring theory and 
observation into agreement within 2a bounds, while remaining in the linear 
regime. Conversely, in the "GUTS" model an increase of a by about 0.04%, 

1. e. A In a = 4 x 10~^, brings theory and observation into agreement within 
la bounds. Considering the variations of fundamental parameters in the three 
scenarios Eqns. (j40l47l49p . the behaviour of the weak scale ((/>) and fermion 
masses is decisive for the variation of abundances. 

6.5 Nonlinear variation of abundances in GUT scenarios 

The unified models discussed in the previous section suggest that it is possible, 
and may even be natural, to obtain a large negative variation in the ^Li abun- 
dance, and considerably smaller variations in other measurable abundances: 
positive in the case of deuterium and negative for ^He. Agreement between 
theory and data in all three abundances could then be possible for a narrow 
range of values in the variation of fundamental parameters, and such scenarios 
could be tested by more accurate abundance measurements. However, the re- 
quired fractional variation in ''Li is so large (a factor two or more in lyLi) that 
a linear analysis using matrix multiplication may be inaccurate. 

We may improve the analysis in specific cases by including the nonlinear 
relations between nuclear parameters and abundances. This is implemented 
simply by running the numerical integration code with the appropriate values 
of nuclear parameters, where the dependence on nuclear parameters was de- 
tailed in Section 13.11 This method would be impractical to investigating the 
full parameter space: nuclear parameters span a 13-dimensional space (or 12- 
dimensional if ^Li is neglected). It is only practicable if the dimensionality of 
the parameter space is reduced by applying unification relations. In principle 
we could also attempt to estimate the nonlinear dependence of nuclear parame- 
ters Xi on fundamental parameters Gk, but this involves additional theoretical 
uncertainty. For the unified models considered here, the fractional variations 
in Xi remain small, well below 0.1. A linear approximation for the relation be- 
tween nuclear and fundamental parameters is therefore appropriate. The main 
nonlinear effects enter at the level of nuclear reactions. 

The nuclear parameters affecting most the large variation in ^Li abundance 
are mainly the deuterium and ''Be binding energies, with the ^He and ^He bind- 
ing energies playing a smaller role. A decrease of Bd causes BBN to happen 
later, which means that the nucleon density is lower and reaction rates smaller. 



22 



ATi/r| 



-0.066 0.000 0.066 




-0.0004 -0.0002 0.0000 0.0002 0.0004 

Aa/a 



Figure 1: Variation of primordial abundances with a in three GUT scenarios. 
Green hues (with smallest variations of abundances) show the first "GUTl"; 
red lines (with large positive slope of Y^u) the second "GUT2"; and blue lines 
(with large negative slope of l^Li) the third "GUTS". Highlighted regions give 
the observational la limits. Error bars indicate the standard BBN abundances 
with theoretical la error |70j, for three different values of r] about the WMAP 
central value, as indicated on the upper horizonal axis. 
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The abundances of ^ > 4 elements are rate-limited and thus decrease with 
decreasing Bj). This accounts for about two-thirds of the change in l/Li- In 
addition, the cross-section of the ^He(a,7)^Be reaction depends strongly on the 
Q-value, hence on the ^Be binding energy. Both these effects are computation- 
ally under control, therefore we believe that the specific nonlinear dependence 
in the scenarios we consider is well estimated within our code. 

We show in Figure [2] the primordial abundances including nonlinear effects, 
i.e. without using a linear approximation for the relation between Ya and Xj. 
For our three GUT models we find a slightly different behaviour of the ^Li 
abundance, which now has an approximately power-law dependence on varia- 
tion of a. It is only slightly more difficult to bring the present observational 
abundances into agreement with standard BBN and the WMAP determination 
of 7] ; still, if we allow a variation of 0.00045 < Alna < 0.0005 in the "GUTS" 
model, the predicted abundances are all very close to the la allowed regions. 

7 Conclusions 

We have developed a systematic method to relate cosmological variations in the 
underlying parameters of particle physics to the primordial isotope abundances 
produced by BBN. The main advantage of the method is that we are able to 
vary every parameter independently, both at the level of fundamental Standard 
Model parameters and of nuclear physics parameters, thus we are not dependent 
on any particular theoretical model which enforces particular relations between 
the variations. 

The method proceeds by defining two linear response matrices. The first, C, 
encodes the change in predicted abundances produced by small variations away 
from the current values of nuclear physics parameters which enter the BBN 
integration code. These parameters (in units where the nonperturbative QCD 
scale is constant) comprise the gravitational constant, fine structure constant, 
neutron lifetime, electron, proton and neutron masses, and binding energies of 
A < 7 nuclei. The dependences of nuclear reaction rates on these parameters 
are also implemented insofar as they are calculated within some effective theory. 
One notable result is that the ''Li abundance depends heavily on the binding 
energies of ^He, ^He and ''Be. 

We also investigated possible further effects of variations in nuclear reaction 
rates on predicted abundances by varying each rate (i. e. thermal integrated 
cross-section (crv)) separately by a temperature-independent factor. We find 
that the ^He abundance is insensitive to nuclear rates, and only eight reactions 
could lead to significant variation of the D, ^He or ''Li abundances. Also in many 
cases the dependences on (av) are small and probably subleading compared to 
other known effects of varying nuclear parameters. 

The second matrix F relates variations in nuclear parameters to the fun- 
damental parameters of particle physics, comprising the gravitational constant, 
fine structure constant, Higgs vacuum expectation value, electron mass, and 
the light (up and down) quark masses. At this point theoretical uncertainty 
enters into the relation between quark masses and nuclear binding energies. We 
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Figure 2: Variation of primordial abundances with a in three GUT scenarios 
including nonlinear effects. Labels as in Fig. 1. 
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parameterise the dependence of binding energies on the pion mass (and hence 
on hght quark masses) by the deuteron binding, which has been treated by a 
systematic expansion in effective field theory. 

The resulting fundamental response matrix R = CF allows us, first, to 
bound the variations of the six fundamental couplings individually, some bounds 
being at the percent level. We then define a "sensitivity matrix" that shows 
which observational data give the best determination of each fundamental pa- 
rameter, and demonstrate the use of the inverse sensitivity matrix to find the 
variations in fundamental couplings required by any given change in primordial 
abundances. 

We can also bound correlated variations affecting many couplings at once: 
we consider three simple scenarios motivated by grand unification of gauge 
couplings. Of these, one allows us to fit observed D, ^He and '^Li abundances 
within 2a bounds, given a variation /S.a/a = —2 x 10~^ away from the present 
value; another fits these observational abundances within Icr bounds, given a 
variation Aa/a = 4 x 10"'^. 

Progress in the field requires both observational and theoretical improve- 
ments. Both statistical and systematic errors in abundance measurements could 
be improved, for example observations to better determine the nature of sys- 
tems where ^He is measured [_25j, or stellar modelling to test possible solutions 
of the ^Li problem. On the theoretical side the relation between quark masses 
and nuclear physics remains unclear beyond the level of the two-nucleon sys- 
tem: the largest uncertainty in our BBN bounds arises from the poorly known 
dependence of the binding energies on the fundamental couplings. 

BBN is already the most powerful probe of fundamental "constants" in the 
early Universe, and precision bounds may well be obtained, given continued 
efforts in observation and theory, to rule out or confirm the presence of a cos- 
mological variation. 
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A Implementation of the BBN code and reaction 
rates 

A.l Units, changes to the code and general remarks 

The code is written in terms of units MeV and 10^ Kelvin. We have to make 
a choice of how to define these units, in a context where ratios of dimensionful 
quantities may be varying. We choose to define units such that the QCD scale 
Ac is constant: thus Ac is always the same number of MeV, and this energy 
scale always corresponds to the same temperature in Kelvin. This simplifies the 
treatment of nuclear physics and QCD; we discuss below in Section fA.SI what 
must be done to translate to other units, or to a unit-free formulation. All 



26 



numerical constants and conversion factors in the code have been replaced with 
variables which are given their values globally, allowing a consistent variation 
to be implemented. 

For the n <^ p reaction rates, instead of analytic approximations [15] we 
implement exact formulae [72j incorporating also zero-temperature and ther- 
mal Coulomb corrections [73]. The default baryon density is taken from the 
WMAP3 determination ry = 6.1 x 10"^'^ ^74j. Note that this determination of 
rj assumes that the values of couplings and mass scales are the same as today. 
Since we keep as a parameter, the possible effect of varying couplings at the 
time of CMB decoupling could be incorporated into our analysis. 

Current measurements of the neutron lifetime are strongly inconsistent, the 
most recent being 878.5 it 0.7 s |75j compared to the 2006 PDG world average 
885.7 lb 0.8s [76|. Since this discrepancy is less than 1% and no abundance 
depends very strongly on r„, we take the lifetime to be the PDG value for the 
purpose of calculating the leading dependence of abundances. 

Starting with the Kawano 1992 code [H] , we replaced all (rounded) numer- 
ical constants in the code by functions of natural "constants", which are glob- 
ally defined as variables. Thus, variations of constants are treated consistently 
within the code. Thanks to increased computational power we can remove all 
significant sources of numerical inaccuracy: these were for example 32 bit in- 
ternal precision of floating point numbers, simplified integration routines, large 
time steps to name just a few. Better numerical accuracy allows us to study the 
behaviour of the abundances under very small changes of "constants" , which is 
essential to derive the linear variation about the standard prediction, and thus 
the derivatives of abundances, without numerical ambiguity. 

A. 2 Fits for charged particle reactions 

Reaction rates for charged particles (with atomic numbers Zi^2 in the initial 
state) arise from a thermal average of a cross-section which in the absence of 
resonances is the product of the Gamow factor and an "S-factor" : 

-27rj? 

a{E) = S{E)-— (51) 



where fj = aZ\Z2\J ^ijlE and is the reduced mass. The S-factor may be 
expanded in a Maclaurin series to quadratic order in energy, which is usually 
sufficient to account for any smoothly- varying dependence. However some cross- 
sections are fit with an additional exponential term 5(0)6"^^ |77j . In addition, 
non-resonant terms may be multiplied by a cutoff factor /cut = e"^"^/"^''"'^ , 
where Tcut has been argued to be proportional to [30l [77]. 

Where the cross-section as a function of energy shows one or more reso- 
nances, they contribute to the thermal averaged rate as 

{av),,, = 5(r)e-^/^ (52) 

where g{T) and E are fitting parameters corresponding to the shape and posi- 
tion of the resonance. Usually a power-law is taken for g{T), thus g{T) = dE^. 
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In principle one should consider the variation of the resonance parameters if this 
term is significant. But since the major contribution to the resonance energy 
E probably arises from A^, the invariant scale of strong interactions, which we 
take as our (non- varying) unit, it seems a reasonable first guess to keep the 
resonance parameters fixed. 

In order to treat the a-dependence of reaction rates consistently, so far as 
it can be calculated, we must implement it at the level of the cross-section a. 
Hence the NACRE formulae [17j fitted at the level of the thermal averaged 
cross-sections {av) are not suitable. Instead we used the functional forms of 
rates from [4:8\ H0| . which correspond to a definite S-factor expansion for a, 
in addition to resonant terms. The free parameters, primarily the expansion 
coefficients of S{E) and resonance parameters, were then fit, in most cases to 
reproduce the NETGEN rates [68] as closely as possible. We also checked that 
the resulting cross-sections are consistent with experimental data. In the case 
of c?(a, 7)^Li we found a set of parameters which seems to fit the experimental 
cross-section at low energies [l9] better than the NACRE fit. But note that 
this cross-section is not measured directly, rather it is derived from experimental 
data under various assumptions, which should be more carefully investigated 

m- 

Replacing the NETGEN rates in the code with our fitted rates, we obtain 
abundances which are changed as follows: Id differs by —0.3%, YsHe by +0.9%, 
l4He by less than 0.1%, Y/u by +3%. Hence we do not consider this refitting as 
significant, except in the case of the d{a, ^)^Li reaction. Depending on whether 
this reaction was fit to NETGEN, or to the cross-section values of [49j, we 
found a ®Li abundance larger by a factor of 1.02, or 3.3, respectively. Given 
the unclear observational status of ^Li this discrepancy is not currently worth 
pursuing. 

A. 3 Scaling of dimensionful nuclear parameters 

The variation of nuclear binding energies and reactions involve (in most cases) 
only two mass scales: the QCD invariant scale Ac, and the average light quark 
mass rh, which indicates the departure from the chiral limit rriq 0. In the chi- 
ral limit the dependence of binding energies and strong interaction cross-sections 
becomes extremely simple: all dimensionful parameters are simply proportional 
to a power of Ac. Switching on the quark masses, one obtains a finite range 
for pion-mediated interactions, which may greatly affect static and dynamical 
properties of nuclei. Also, the masses of all hadrons are affected at some order 
in chiral perturbation theory [58j. However, if both Ac and ruq are varied by 
a common factor, while all dimensionless couplings are held constant, dimen- 
sionful quantities involving strong interactions (and to a good approximation 
electromagnetic interactions) scale with some power of this common factor. 
Such a variation of dimensionful parameters is equivalent to a redefinition or 
variation of the units of mass or energy. 

Since we fix our unit of mass to be a constant times Ac, the correct behaviour 
of dimensionful quantities associated with QCD or the strong nuclear force is 
automatic. So long as quark masses are proportional to Ac, such quantities are 
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unchanged. The effects of other dimensionful parameters (0) and Gn are then 
transparent. 

If, however, is formally allowed to vary, the scaling properties of QCD, 
and hence of the strong nuclear force, provide a simple check on the dependence 
of physical quantities on dimensionful parameters. For example, we will write 
the dependence of the deuteron binding energy on the pion mass as 

A In 5d = r A In (53) 

in units where Ac is constant. We may thus rewrite variations of dimensionful 
quantities as, for instance, Aln(i?D/Ac) and obtain 

Aln^D = rAlumTT + (1 - r)AlnAc, (54) 

allowing us to compare results obtained with different choices of unit. For this 
type of equation relating fractional variations of dimensionful quantities, there 
is a simple check on whether it behaves correctly under redefinition of units. 
Each term is a numerical coefficient rj times the fractional variation A In Qj of 
a quantity Qi with mass-energy dimension Di: then the sum of rjDj on each 
side of the equation must match. 



B Observational situation and uncertainties 

One of the biggest success of standard BBN is the matching of theoretically 
predicted and observed primordial abundances for major elements. For a re- 
view of the theoretical and observational status and obstacles see [25]. The 
highest precision measurement is that of the ^He abundance (conventionally 
written Ip); however the actual precision of the determinations and possible 
systematic errors are currently debated [25]. The situation is exemplified by re- 
cent contradictory observational determinations [80j. It was argued in [69j that 
given a range of systematic effects the observational data indicated a primordial 
abundance of 

Yp = 0.249 ± 0.009. (55) 

which we take to be a la range. However, given the probable dominance of 
systematic effects, instead of using 2-cj bounds to determine the range of allowed 
variations, we rather use the "conservative allowable range" of Yp given in [69] 
as 

0.232 <Yp< 0.258. (56) 

The determination of the primordial deuterium abundance follows from a 
small number of observed systems. Recent determinations |8H I82j yield a value 
of 

D/H = (2.8±0.4) X 10"^ (57) 

where the large scatter between values determined from different systems should 
be noted. 

Determinations of the ^He abundance typically have a large scatter. Consid- 
ering also the complex post-BBN development of this isotope, ^He abundance 
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determinations cannot be considered as good tracers for the primordial abun- 
dance |i83j. 

Observations of the ''Li abundance show a plateau at old, metal-poor halo 
stars, suggesting that the plateau value is closely related to the primordial one. 
The most recent determinations of the abundance are quite small: ^Li/H = 
(1.3 ± 0.3) • 10~^° [H] (see also It has been suggested that there are 

unresolved systematic errors relating to the effective temperature of the stars 
[85] which may imply a value as large as ^Li/H = (1.64ib0.3) • lO"^'^. To account 
for this possible systematic, we adopt a value 

^Li/H = (1.5±0.5) X 10-^°. (58) 

Thus the observed ^Li abundance is about a factor of 3 smaller than the stan- 
dard theoretical prediction. 

A possible detection of ^Li was discussed in [78], though not at high sta- 
tistical significance. If the detection is correct, the ^Li abundance is about a 
factor 1000 larger than the SBBN prediction. Given the unclear observational 
status and post-BBN history of the isotope, we do not include ^Li in the final 
analysis. 

Theoretically predicted primordial abundances also come with an error, 
mainly due to cross-section uncertainties. Our numerical procedures do not 
provide error estimates, so we adopt the la ranges from f7U\, using a baryon 
density Q^/i^ = 0.0224 

D/H = (2.61 ±0.04) X 10^^ 
^He/H = (1.03 ± 0.03) x 10"^ 

Yp = 0.2478 ± 0.0002 
"^Li/H = (4.5±0.4) X 10~^°. (59) 
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